Ensemble Averaged Conductance Fluctuations in Anderson Localized Systems 
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We demonstrate the presence of energy dependent fluctuations in the locahzation length, which 
depend on the disorder distribution. These fluctuations lead to Ensemble Averaged Conductance 
Fluctuations (EACF) and are enhanced by large disorder. For the binary distribution the fluctua- 
tions are strongly enhanced in comparison to the Gaussian and uniform distributions. These results 
have important implications on ensemble averaged quantities, such as the transmission through 
quantum wires, where fluctuations can subsist to very high temperatures. For the non-fluctuating 
part of the localization length in one dimension we obtained an improved analytical expression valid 
for all disorder strengths by averaging the probability density. 



In quantum coherent conductors, disorder induces dra- 
matic effects on the conductance. For instance, if the lo- 
calization length, Lc, which is the decay length of the 
wavefunction, exceeds the system size L, the conduc- 
tance vanishes exponentially (Anderson localization) [l| . 
On the other hand, if the system size is smaller than 
Lc, then the conductance exhibits fluctuations in energy 
which take on universal amplitudes and are referred to 
as Universal conductance fluctuations (UCF) \T^ . In one 
dimension, the effect of disorder is particularly dramatic 
because the presence of any amount of uncorrelated disor- 
der will lead to a finite localization length at all energies. 
The localization length, which is an ensemble averaged 
quantity, has typically no fluctuations as a function of 
energy. This is in stark contrast to the conductance for a 
single disorder configuration, where the fiuctuations are 
maximized due to UCF. 

In many quantum wires, such as carbon nanotubes, 
the conductance can be measured as a function of the 
gate voltage, effectively changing the Fermi energy of 
the carriers. As a function of energy, the conductance 
often shows strong fluctuations associated with disorder. 
The amplitude of these fluctuations decays with temper- 
ature [3, kl . This can be understood in terms of ensem- 
ble averaging, since an increase in temperature leads to 
a decrease of the quantum coherence length Iff, due to 
increased inelastic scattering. This in turn leads to the 
ensemble averaging over the quantum size /^ averaged 
L/l^ times for a sample of size L. it is therefore surpris- 
ing that experimentally, the conductance fluctuations of 
carbon nanotubes remain substantial even at high tem- 
peratures (equivalent to ensemble averaging) [3|, IJ] ■ 

In order to better understand possible fluctuations 
even after ensemble averaging, we studied in detail the 
energy dependence of the ensemble averaged quantity Lc- 
Quite strikingly, we observe fluctuations of Lc as a func- 
tion of energy, which depend on the distribution of the 
disorder potential. This is illustrated in Figures [T] and 
[21 where we show the relative variation of the inverse lo- 
calization length (Lyapounov exponent, A) as a function 
of energy for different disorder distributions (Gaussian, 
uniform and binary). For the calculations we used the 



one-dimensional Anderson model, described in more de- 
tail below. 
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FIG. 1: The relative deviation of the Lyapounov exponent for 
the binary and skewed binary distribution. The skewed binary 
distribution function of average zero is given by P{V) = (1 — 
a)5{V - aW) + aS{V + {1 ~ a)W), where we have a = 1/2 
for the binary case (As) and we chose a — 1/3 for the skewed 
case (Ass). The variance a^ is varied between .05 and 1.3 
and depends on W. Xp is the Lyapounov exponent obtained 
from the average probability density. 



The variations (or fluctuations) of the Lyapounov ex- 
ponent are quite substantial for the binary distribution 
and reach up to 80%. The fluctuations are shown with 
respect to the Lyapounov exponent (A^) obtained from 
the average probability density, valid for all disorder 
strengths and derived below. The symmetry of the fluc- 
tuations depends on the symmetry of the distribution. 
Indeed, for a = 1/2 the distribution is symmetric in en- 
ergy, which leads to symmetric fluctuations as seen in 
Figure [1] in stark contrast to the skewed case a / 1/2, 
where the distribution is asymmetric. 

In Figure [2] we show the relative deviations of the 
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FIG. 2: The relative deviation of the Lyapounov exponent 
for the Gaussian distribution (Ac) of zero average and of 
standard deviation a and the uniform distribution {\u with 



-W/2 <V < W/2 and (j'^ = W^/12). 
0.05 and 1.3. 



a is varied between 



Lyapouvov exponent for the Gaussian and uniform dis- 
tributions. Here the fluctuations are strongly suppressed 
and the deviations are of the order of 10% and are very 
smooth when compared to the binary distributions for 
the same standard deviation range. However, the ob- 
served strong fluctuations for the binary distributions, 
lead us to make a comparison with mesoscopic conduc- 
tance fluctuations. 

In general, conductance fluctuations are expressed as 
{SG^) = {{G — (G))^), where (•) denotes the ensemble 
average, i.e., the average over a given disorder distribu- 
tion. The standard universal conductance fluctuations 
(UCF) result gives ^/{SCP) ~ Q.73e'^/h ^] for a quasi 
ID system. In the localized regime these fluctuations 
are suppressed. In order to characterize the fluctuations 
in the Lyabounov exponent seen in Figure [1] and [2l we 
define 
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[X{E) - X,{E)]ME, (1) 



where Xp is the Lyapounov exponent obtained from the 
average probability density with the same standard de- 
viation. The bandwidth {Emax — Emin) is fixed to the 
non-disordered one. 

For a strongly localized system of size L, the ensemble 
averaged conductance is given by (G) ^ exp"^'*'-'". Hence, 
the fluctuations in A lead to fluctuations in the ensemble- 



averaged conductance 5{G) ^ 2{G)L5X, which we coin 
ensemble averaged conductance fluctuations (EACF). We 
studied EACF for different disorder strengths. The re- 
sults are shown in Figm'e[3]for different distributions. 
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FIG. 3: S\^ is shown for different distributions as a function of 
the disorder strength, which is characterized by the variance 
a^. For the skewed binary distribution we used a = 1/3 (the 
same as in Figure [TJ. These fluctuations represent ensemble 
averaged conductance fluctuations (EACF). 



We now describe the derivation of the Lyapounov ex- 
ponent (Xp) for the average probability density, which is 
used in the calculations of the relative fluctuations shown 
in Figure [21 This is a novel analytical approach in order 
to obtain an expression for X{E) valid for all disorder 
strengths. 

The one-dimensional Anderson model 1] is a tight 
binding equation with random on-site potentials Vn and 
given by 



i*„+i-fi*„_i = (^-F„)*, 



(2) 



Here t is the hopping term, which we set to one. In this 
model, an alloy of two elements would be described with 
V„ taken from a binary distribution like in Figure [TJ In 
contrast, if the disorder is due to the surface roughness 
of the substrate, like in the case of a carbon nanotube on 
a silicon oxide substrate, Vn would be given by a more 
continuous distribution, like the Gaussian or uniform one. 
In order to obtain X{E) the main idea is to obtain and 
iterative equation for the probability density and then 
to average it. Hence, assuming real potentials Vn it is 
possible to rewrite ([2]) as 



Pn+1 



E-Vn 



[{E - K 



{E-V^.i} 



1] P. 



{E - Vnf - ,^ ,;" T Pn + [1 - (S - Vn){E - K-l)] Pn-1 + 

\-2{E- Vn)iE - Vn-l) + {E- Vn-lf + 

+ [{E - VnKE - K_2) + 1 - (^ - Vn-l)iE - Vn-2)] Pn-2 + 

I 



E-Vn 



E-Vn-1 

Vn - K^l 

{E-Vn-2} 

Vn-1 - V„ 

Pn~3, 



E - Vn- 



Pn-2 
Pn-l 



(3) 
(4) 



where p„ = ipn'4'n is the probabihty density. Equation 
(|4]) is obtained by using an additional iteration and wih 
be important when considering the average. Interest- 
ingly, this exact (before averaging) iterative expression 
for the probability density depends explicitly on neigh- 
boring potentials, which illustrates the importance of the 
assumption of an uncorrelated disorder potential when 
taking the disorder average. Indeed, the presence of lo- 
cal correlations can lead to an infinite localization length 
in one and two dimensional disordered systems (5| . In our 
method of considering the probability density, p„ instead 
of ipn , phase correlations are not averaged out when the 
disorder average is performed. 

The average of ^ can be taken trivially by assuming 
uncorrelated disorder with {VnVm) = o'^Sn.m and noting 
that {piVj) = {pi){Vj) for i > j, since p„ does not depend 
on Vn when pn is obtained iteratively using initial con- 
ditions for n = 0, 1,2,3. The term in curly brackets of 
Pn in © is the only product which cannot be separated 
from Pn when averaged, which is the reason we had to 
iterate this equation one more time in order to obtain 
equation ([4]), where the coefficient in front of the curly 
bracket term now averages to zero. 

For a potential of average zero we thus obtain the fol- 
lowing iterative solution for the average probability den- 
sity: 



{pn+l) = {E''+a^-l){pn) + {l-E''+<T''){pn-l) + {pn-2)- 

(5) 
The leading dependence of (p„) can now be extracted 
by evaluating the eigenvalues {Ci, $2, Ca} of the character- 
istic transfer matrix determined by ([5]). The correspond- 
ing Lyapounov exponents are 



A* = 2l^°s(l6 



(6) 



where the factor 1/2 comes from the definition of the Lya- 
pounov exponent in terms of the probability density in- 
stead of the wavefunction. For cr = all three Lyapounov 
exponents are zero, which is characteristic of the plane 
wave solution for an energy within the band —2 < E < 2. 
Defining Ai < A2 < A3 we have for cr > that Ai > 0, 
which implies that all states are localized. An exception 
occurs at ii^ = 0, where Ai = and A2 > 0. However, 



the solution corresponding to Ai violates the positivity 
requirement of p„. Hence, for energies close to i? = we 
arbitrarily define Xp as the average of the two smallest 
Lyapouov exponents, instead of just the smallest expo- 
nent, i.e., Xp = (Ai + A2)/2. Away from the band center, 
corresponding to \E\ > o''^/4, we have Ap = Ai = A2. 
Here the smallest Lyapounov exponent Ai correctly de- 
scribes the dependence. In the limit of small disorder 
{a < 1) this leads to Ap ~ (l/2)crV(4 - E^) = XriE), 
where At corresponds to the standard result by Thouless 
[6|, which was obtained by averaging the Green's func- 
tion to second order. For |£;| > 2 -I- 3a'^/^/A {a < 1) 
we have Ai y^ A2 and the solution corresponding to Ai 
violates the positivity requirement of p„, which is why 
we define Ap = A2 in that range, which is outside of the 
non-disordered band. The behavior of the Lyapounov ex- 
ponents corresponding to the three eigenvalues are shown 
in the inset of Figure [Hand are similar to the ones studied 
for continuous disordered systems (7[. 

For the special energy at the band center, the correc- 
tion for the small disorder expansion [sj, |9| is given by 
X'^-^iE = 0) ~ 0.91A'^(0), where A^"^ is the correc- 
tion due to the Kappus-Wegner anomaly. This anomaly 
at the band center can also be seen in the context of 
our approach using the probability density, since expres- 
sions ([3] and [4|) are ill defined when taking the average 
at ii^ = 0. At the band edge the first order correction 
to the Thouless result was obtained recently by using a 
classical two-dimensional Hamiltonian map [lOj . 

The anomalies at the band center and at the band 
edges can be easily identified in Figure [2] for the Gaus- 
sian distribution, where they appear as bumps for small 
disorder. Most previous analytical approaches are based 
on small disorder expansions and cannot be applied to 
high disorder. This is in contrast to our probability den- 
sity approach, where we didn't assume small disorder, 
and we indeed obtain the correct large disorder limit 
Ap = A2 — A3 ~ log(cr) with a Lyapounov exponent 
independent of E 
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Numerically, the Lyapounov exponent was obtained by 
evaluating the eigenvalues of the product of transfer ma- 
trices obtained from equation ^, 
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(7) 



and then taking the self-averaging hniit where we used 
TV > 10^ 

A- lim llog(max{|^f|}) = - hm 1 log(min{|ef |}). 

(8) 
The quaHty of our expression for Xp can be seen in 
Figure m where we show the relative deviation to the nu- 
merically obtained Lyapounov exponents Xg,u- The de- 
viations are of the same order as the dependence on the 
distribution function (Gaussian versus uniform). This 
shows that any further improvement of the analytical 
expression for X{E) needs to depend on the distribution 
function explicitly (for example by including higher mo- 
ments) . For the special case of the Cauchy distribution, 
where the second moment does not exist, an exact an- 
alytical expression for the Lyapounov exponent can be 
found [i2. 
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FIG. 4: Lyapounov exponents {Xq, Xp and At) as a function 
of energy for different disorder strengths. In the inset Ai, A2, 
and A3 are shown together with Xq and Xp. 



We also compared our expression Xp in Figure [4] for 
a large range of disorder strengths and find that the fit 
is less than a factor two off for all energies and disorder 
strengths, in contrast to the Thouless expression (At) 
also shown, which deviates substantially at large disor- 
der. 

The ensemble averaged fluctuations, which is the main 
result presented here, were obtained for one of the sim- 
plest models showing Anderson localization, namely the 
random one-dimensional tight binding equation. How- 
ever, Anderson localization is very general and can be 
observed, among others, in the propagation of light in dis- 
ordered media, in phonon and plasmon modes, in quan- 
tum chaotic systems ^IBj . in Bose Einstein condensates 



1J|, and even in neutron propagation [15|. We there- 
fore expect that similar EACF also exist in these sys- 
tems, since the equations describing these are very sim- 
ilar. For dimensions higher than one, only very few an- 
alytical expressions exist, including a recent result valid 
for small disorder, which has been derived for quasi one- 
dimensional systems [16| and that can be expressed as 
a convolution of the one-dimensional case. We there- 
fore expect that at higher disorder strengths, distribu- 
tion dependent fluctuations will also appear in quasi one- 
dimensional systems, which are relevant for the experi- 
ments on quantum wires. More striking examples might 
be systems composed of alloys, since the potentials of 
a two component alloy would be described by a binary 
disorder distribution, which shows the strongest EACF. 
Indeed, in InGaAs alloys, conductance fluctuations in a 
quasi one-dimensional geometry were observed at high 
temperatures [17|. However, the study of EACF in higher 
dimensions is beyond the scope of this work. 
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